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The first resolved imaging of milliarcsecond-scale jets in 
Circinus X-l 
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ABSTRACT 

We present the first resolved imaging of the milliarcsecond-scale jets in the neutron 
star X-ray binary Circinus X-l, made using the Australian Long Baseline Array. The 
angular extent of the resolved jets is ~ 20 milliarcseconds, corresponding to a physical 
scale of ~ 150 au at the assumed distance of 7.8 kpc. The jet position angle is relatively 
consistent with previous arcsecond-scale imaging with the Australia Telescope Com- 
pact Array. The radio emission is symmetric about the peak, and is unresolved along 
the minor axis, constraining the opening angle to be < 20°. We observe evidence for 
outward motion of the components between the two halves of the observation. Con- 
straints on the proper motion of the radio-emitting components suggest that they are 
only mildly relativistic, although we cannot definitively rule out the presence of the 
unseen, ultra-relativistic (r > 15) flow previously inferred to exist in this system. 
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1 INTRODUCTION 

Circinus X-l (Cir X-l) is one of the few confirmed Galac- 
tic neutron star X-ray binaries with resolv ed radio jets 
which have been studied from sub-arcsecond (ICalvelo et al 



20 ll[ ) through arcsecond (|Fender et al] 19981: Tudose et a! 
20081 ) and out to arcminute scales (IStewart et all 1 19931 : 



Tudose et al.ll2006l ) , where they have inflated a synchrotron- 
emitting radio nebula. Cir X-l is also the only confirmed 
neutron star system with resolved X-ray jets, or iented along 
the s ame position angl e as the radio jets |Heinz et al.1 
120071 : ISoleri et al.1 l2009ri. The system comprises a neu- 
tron star primary (jLinares et"all l20ld ) accreting from a 
less-evolved companion whose nature is relatively poorly 
constrained. The binary has an eccentric 16.6 day orbit 
l|Kaluzienski et all 19761), with enhanced accretion close to 
periastron (jMurdin et al ] ll98d : lHavnes et aflll98d ) leading 
to flaring events in the radio, infrared and X-ray bands. The 
distance to Cir X-l is not well-constrai ned. From the peak 
luminosity of its Type I X-ray bursts, Ijonker fc Nelemansl 
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l|2004 ) estimated a distance in the range 7.8-10.5 kpc, in 
agreement wit h the earlier Hi a b sorpt ion distanc e of > 8 kpc 
determ ined bv lGoss fc Meboldl (|l977l ). However. Ilaria et al.1 
|2005l ) claimed a lower value of 4.1 kpc from model- fitting 
of X-ray spectra that suggested a low hydrogen column to- 
wards the source. 

The orientation of th e system remains controversial. 
The presence of X-ray dips (Shircy ct al. 1999) and the spec- 
tral changes on egress (|Brandt et al.1 1996J). as well as the 
P-Cy gni profiles of detected X-ray lines l|Brandt fc Scmibd 
l200d ). argue for an edge-on accretion disc. However, the ex- 
tremely high borentz factor (r > 15) inferred from the time 
delay between a flare in the radio core and subsequent down- 
stream flaring in the jets would require the radio jets to b e 
inclined very close to the line of sight (|Fender et al.|[2004l ). 
Reconciling all these constraints would require a bent jet or 
a severe misalignment between the jet a xis and the orbit al 
plane (as discussed for other systems by lMaccaronejl2~002t ) . 

The ultra-relativistic flow inferred by iFender et al.1 

|200ll ) depends on the assumption of a causal relation be- 
tween flaring events in the core and those subsequently ob- 
served downstream in the jets. However, the sparse time 
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coverage of the existing radio observations has thus far pre- 
cluded independent confirmation of this result. With the 
high angular resolution available using very long baseline 
interferometry (VLBI), we can directly measure the proper 
motions of individual jet components, hence constraining the 
jet Lorentz factor and inclination angle. However, compact 
radio emission is seen only during the radio flares immedi- 
ately following periastron (|Moin et alj|201ll ), and the low 
surface brightness sensitivity of VLBI arrays constrains the 
observations to periods of strong radio flaring in the system. 

From 1975-1985, the radi o flares regular l y rea ched lev- 
els in excess of 1 Jy, allowing IPreston et all (l983h to con- 
strain the angular size of the flaring component to lie be- 
tween 1.5 and 15 mas at 2.3 GHz. After 1985, the source 
entered a more quiescent phase with flares peaking at mjy 
levels, and further activity was not seen again until 2006, 
when seven flares ranging from 0.2-1 Jy were detected at the 
Hartebeesthoe k Radio Astronomy Observatory (HartRAO; 
lNicolsonll2007l ). This triggered new VLBI observations of Cir 
X-l, and in the first Southern H emisphere electronic VLBI 
experiment, IPhillips et all (|2007n detected a single, scatter- 
broadened component of size 60 ± 15 mas at 1.6 GHz. The 
source was also unresolved in the 1.4- and 1.7- GHz e-VLBI 
observations of iMoin et al] l|201ll ). 

In 20 10 May, enhanced X-ra y activity was detected from 
Cir X-l (jNakaiima et al] 120101 ). Follow-up radio observa- 
tions were carried out at HartRAO (Nicolson, priv. comm.) 
and subseq uently at the Austr alia Telescope Compact Ar- 
ray fATCA; ICalvelo et al.|[201ol ). measuring flux densities of 
order several hundred mjy, an order of magnitude higher 
than seen at the ATCA over the past 15 y. This prompted 
us to trigger new VLBI observations at a higher frequency 
and hence improved angular resolution, in an attempt to de- 
termine the milliarcsecond-scale morphology of the source. 



2 OBSERVATIONS AND DATA REDUCTION 

We conducted a target-of-opportunity observation of Cir X- 
1 using the Australian Long Baseline Array (LBA) from 
02:21-16:41 UT on 2010 July 28 (orbital phase 0.046-0.082 
according to the latest ephemeris from G.D. Nicolson, priv. 
comm.). We observed at the relatively high frequency of 
8.4 GHz, both reducing the angular broadening caused by 
interstellar scattering and increasing the achievable angular 
resolution. The array comprised the ATCA (in its hybrid 
H168 configuration), Ceduna, Hobart, Mopra and Tidbin- 
billa telescopes. The ATCA was used in tied-array mode, 
phasing together the four antennas CA01, CA02, CA03 and 
CA05 (maximum baseline of 192 m). No fringes were de- 
tected from Tidbinbilla, reducing our effective LBA array to 
four stations. The resulting aii-coverage is shown in Fig. [1] 
The observations were conducted in dual-polarization for- 
mat with 64 MHz of bandwidth per polarization. 

Th e data were correla t ed us ing the DiFX software cor- 
relator dDeller et al] 120071 , | 20 111 ), and reduced using aips 
i|Greisenll2005M ~ The source amplitudes were determined us- 
ing the phase calibrator source J1515-5559, assuming an ap- 
proximate flux density of 1.7 Jy from recent ATCA monitor- 
ing (with an estimated error of ±20 per cent), an estimate 
confirmed to within the uncertainties from the ATCA data 
taken in VLBI mode during the LBA observations. Since 
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Figure 1. uw-coverage of our observations. Red points indicate 
the first half of the observations, and black points the second half. 



Cir X-l was sufficiently bright, fringe fitting was performed 
on the target source. I maging and self-calibration were the n 
carried out in Difmap (|Shepherd. Pearson fc Tavlorll 19941 ). 



3 RESULTS 

Cir X-l was resolved by the LBA into a symmetric, jet-like 
structure that evolved over the course of the observations 
(Fig. 0. Although our minimum and maximum baselines 
were 113 and 1702 km, respectively, the failure at Tidbinbilla 
deprived us of intermediate baselines of order a few hundred 
kilometres. This created a large hole in the uv-plane (Fig.[T]), 
such that we were only sensitive to structures on scales of 4- 
9 and 67-250 mas. Our single short baseline (between ATCA 
and Mopra) was consistent with a point source whose flux 
density gradually decayed from 210 to 80mJybeam~ . 

The changing flux density over the 14.3 h observation 
violates one of the fundamental assumptions of aperture 
synthesis. In the absence of morphological variations, this 
would create amplitude errors in the image plane. However, 
since our w^coverage is already so sparse, we cannot com- 
pensate by only imaging short snapshots of data in which 
the amplitudes are relatively stable. Tests indicated that the 
data could not be binned more finely than into two separate 
halves without severely compromising the imaging fidelity. 
A comparison of the data from the two halves of the ob- 
servation (Fig. [2} shows a fading of the core jets, and the 
appearance of two new outer components in the second half. 
If these new components correspond to outward motion of 
the core jets seen in Fig. [2ja) , the combined proper motion 
of the two components is ~ 35 mas d -1 , significantly differ- 
ent from the 400 mas d -1 inferred by iFender et al] (|2004l ) 
from the time delay between core and lobe flaring events. 
The expansion of the jets corresponds to the appearance of 
structures on larger scales, reflected in the shift to shorter 
spacings of the first minimum in the amplitudes of the visi- 
bility data projected along the jet axis (Fig. [3] top) . 
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(a) First half of the observations (rms noise 
1.51 mjy beam -1 ). Total flux density of the jet inside 
the 3c contour is 146 ± 14 mjy. 




(b) Second half of the observations (rms noise 
0.75 mjy beam -1 ). Total flux density of the jet inside 
the 3<t contour is 104 ± 8 mjy. 



Figure 2. LBA images of Cir X-l on 2010 July 28. Contours for both images are at ±3 X (%/2) n times the rms noise level, where 
n = 0,1,2,.... The colour bar shows the radio brightness in units of mjy beam -1 , with the scale being the same for both images. 
Restoring beam size is 3.4 X 2.9mas 2 in P.A. —69.1°. The source is resolved into a continuous jet, which expands from an extent of 
9.3 mas to 19 mas between the two halves of the observing run. 



To evaluate the effects of the changing source flux den- 
sity and morphology combined with our relatively sparse 
wv-coverage, we performed some simulations within AIPS. 
We subtracted the clean components from our best image 
(Fig. [2]) in the ?ti>-plane, and then added various source mod- 
els back in to the uv-data,. A decaying point source at the 
phase centre as seen on the ATCA-Mopra baseline could 
not reproduce the symmetric extended structure along the 
observed position angle. We therefore conclude that the ob- 
served structure is not an artifact caused by amplitude errors 
combined with limited m>coverage. Adding a second, non- 
varying component along the position angle of our observed 
jet showed that we could not reproduce the observed sym- 
metric structure with an asymmetric, one-sided jet model. 
Adding instead a single, ballistically moving jet component 
created a one-sided, smooth jet, with a reduced surface 
brightness owing to the smoothing of the emission along the 
component tr ajectory, and a slig htly curved jet locus (as 
also found bv lTingav et al.lll995l ). Artifacts also appeared 
as stripes parallel to the jet trajectory, reminiscent of the 
noise structure in Fig.[2ja). From these simulations, we con- 
clude that while the relatively continuous appearance of the 
jets could be created by discrete, moving components, the 
symmetric profile we observe is real, as it cannot be repro- 
duced using any one-sided source model. The reality of the 
symmetric structure in the first half of the observations is 
also demonstrated by the flat and zeroed phases of the uv- 
data when projected along the jet axis (Fig. [2]). Some slight 
asymmetry on scales of 7-10 mas is visible in the second half. 

Finally, we tested whether the extra components seen 
in Fig[2jb) could have been detected with the tii>coverage 
available in the first half of the observations. Substituting 
in the uihpl&xie the clean components of Fig. (2jb) for those 
of Fig. [2ja), we found that the extra emitting components, 



if present in the first part of the observations, would have 
been clearly visible in the images. Either the jets have moved 
outwards between the two halves of the observation, or new 
components have become visible, via energisation of the elec- 
trons or amplification of the magnetic field. 

The observed jets are elongated along a position angle 
112° E of N, in good agreement with the position angle of 
the arcseco nd-scale jets detect ed in both the radio and X- 
ray bands. I Fender et alj (| 19981 ) derived a position angle of 
110 ± 10 ° for the extended a rcsecond-scale radio emission, 
although iTudose et al.l (|2008h found this position angle to 
vary between epochs, with a mean of 129 ± 13°, although 
there was no une quivocal evidenc e for precession. The deep- 
est X-ray images (|Sell et al.ll2O10h show unresolved emission 
at a position angle 140° E of N, and extended diffuse emis- 
sion from 95-155° at a distance of 20-50" from the core. 
From the sli ghtly different orient ations seen on different an- 
gular scales. ICalvelo et alj (|201ll ) suggested either a varying 
degree of precession with time or bends in the jet. 

We find no evidence for angular broadenin g of the radio 
emission, as reported bv lPhillips et al] (|2007T ) at the lower 
frequency of 1.6 GHz. Their measured angular size of 60 ± 
15 mas, scaled by the expected A 2 ' 2 dependence corresponds 
to 1.6 ± 0.4 mas at our observing frequency of 8.4 GHz, well 
below our beam size of 3.4 x 2.9 mas 2 . 



4 DISCUSSION 
4.1 Morphology 

The symmetry of the jets is initially reminiscent of the com- 
pact jets observed in the hard state of the black hole sys - 
tem GRS 1915+105 (|Dhawan. Mirabel fc Rodrigue j [200ol ) , 
although we note several caveats to this interpretation. The 
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Figure 3. Visibility amplitude and phase as a function of pro- 
jected ito-distancc along the position angle of the jet axis. Plotted 
points show the data from the second half of the observation. Red 
and black lines show the clean component models for the first 
and second halves, respectively. The total amplitude decreases 
between the first and second halves, and the first minimum shifts 
to shorter projected baselines, implying that the overall structure 
has expanded. The relatively flat phases as a function of projected 
baseline length show the symmetry of the jets about the peak. 

evolution of the source structure over the duration of the 
observing run (Fig. [2j implies that this is not a persistent, 
compact jet. Second, while we have no spectral information 
for this particular flare, previous multi-frequency observa- 
tions of Cir X-l have shown that the spectra evolve from 
optica lly thick to optical l y thin over the cou rse of an out- 
burst (|Fender et al.ll2005l ; lCalvelo et alj|2010h . and are usu- 
ally optically thin by orbital phase 0.05 (G.D. Nicolson, priv. 
comm.). Thus it is likely that the observed emission is caused 
by the outward motion of discrete, expanding, optically-thin 
eject a, as seen during outbursts of black hole X-ray binaries 
(e.g. iMirabel fc Rodriguej[l993 : iTingav et aDll995T l. 

Since the observations were not phase-referenced, we 
have no absolute positional information. Also, the astromet- 
ric parameters of Cir X-l are unknown, so the predicted loca- 
tion of the binary system cannot be accurately determined. 
It is therefore uncertain whether the peak brightness corre- 
sponds to the location of the binary system, or rather to the 
surface of optical depth unity in the approaching jet. The 
latter possibility would require a symmetric profile for the 
approaching jet, and for the receding jet to be hidden either 
by Doppler de-boosting, or by free-free absorption within a 
few hundred au of the central binary. However, the appar- 
ently symmetric nature of the expansion seen in Fig.[2]tends 
to argue against this scenario. The most plausible explana- 
tion for the observed symmetry is therefore an intrinsically 
slow jet such that the Doppler factors of approaching and 
receding jets are comparable. In that case, the jets are only 
mildly relativistic and inclined close to the plane of the sky. 

4.2 Jet speed 

An ultrarelativistic flow with a proper motion of 400 mas d~ 
would be smeared over 68 beams over the course of our 
14.3 h observation. Unless the ejecta were intrinsically ex- 
tremely bright (> 112 mjy), this would reduce the surface 
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Figure 4. Intrinsic jet speed as a function of inclination angle, 
given the proper motion constraints of /^ a pp + Mrcc ~ 16masd — 1 
(black curves) and 35 mas d -1 (red curves). Curves are plotted for 
the assumed distance, as well as for the lower and upper bounds of 
the commonly accepted distance range for the source. Our proper 
motion constraints are both compatible with a mildly relativistic 
jet velocity at any moderate inclination angle, unless the distance 
is at the upper end of the estimated range. 

brightness bel ow our 3cr detection threshold. However, in 
the scenario of I Fender et al] (|2004h , the postulated ultrarel- 
ativistic flow was unseen, but inferred from the brighten- 
ing of radio-emitting components downstream in the jets as 
they were energized by this faster outflow. A similar model 
was used to explain the brightening of radio-emitting lobes 
moving away from the neutron star X-ra y binary Sco X-l 
ijFomalont. Geldzahler fc Bradshawll200lf ). While the lobes 
in Sco X-l moved at speeds of 0.32-0. 57c, energy appeared 
to move from the core to the lobes at > 0.95c. A similar 
phenomeno n has also been seen in the jets of active galactic 
nuclei fe.g. ITingav et al] 1 19981 ). Thus we might not expect 
to directly detect radio emission from any ultrarelativistic 
flow. 

Since compact radio emission fro m Cir X-l is onl y 
present at and shortly after periastron (|Moin et alj|201lh . 
we can assume that the jets were energized by the flare 
at orbital phase zero. If an ultrarelativistic flow causes the 
brightening of the extra components in Fig.EHb), we can use 
the 1.2 d delay between phase zero (from our ephemeris) and 
the time of the observations to infer its speed. The observed 
size of 19 mas then gives a combined proper motion of ap- 
proaching and receding jets of ~ 16 mas d" 1 . This value is 
even smaller than the combined proper motion of 35 mas d 
calculated from the expansion seen in Fig. [2] suggesting that 
the actual jet ejection might have come slightly later than 
phase zero. Although these proper motion estimates do not 
allow us to uniquely determine the jet speed and inclination 
angle, for a given distance we can determine the jet speed as 
a function of inclination (Fig. |4]). Unless either the jets are 
inclined extremely close to the line of sight, or the distance 
is at the upper end of the allowed range and the estimate of 
35 mas d -1 is most representative, the intrinsic speed of the 
radio-emitting components is only mildly relativistic. 

A flow travelling at 400 mas d _1 should by contrast have 
propagated out to 0.5 arcsec by the end of our observations. 
Imaging out to 0.8 arcsec from the field centre showed no ad- 
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ditional emission downstream in the jets, so any downstream 
knots energized by an ultrarelativistic flow must have been 
sufficiently diffuse to be resolved out. Should an ultrarel- 
ativistic flow exi st, the jets must b e aligned within 5° of 
the line of sight (jFender et al.|[20o3 ). While reconciling this 
constraint with the symmetric morphology discussed in Sec- 
tion l4,ll is difficult, the uncertainties in our interpretation are 
sufficient that we cannot definitively rule out this scenario. 
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4.3 Opening angle 

Owing to our sparse u^coverage coupled with intrinsic 
source variability, we cannot derive meaningful constraints 
on the opening angle from the visibility data. In the image 
plane, the jets are not significantly resolved transverse to the 
jet axis. From the length of the jets and the beam size of 

3.4 x 2.9 mas 2 in PA -69.1°, we constrain the opening angle 
of the jets to be < 20° . This allows us to rule out poor colli- 
mation as an explanation for the hal f-opening angles of 35° 
found for the X-ray caps observed bv lSell et all l|2010l ). The 
remaining options put forward by the authors are precessing 
jets or a jet axis aligned very close to the line of sight. The 
variation in the position angle of the extended radio emission 
and the symmetry of the milliarcsecond-scale structure (Sec- 
tion [3]) when compared with pre vious arcsecond-scale radio 
observations (|Calvelo et al . 2011) would seem to support the 
precession scenario, although the standard deviation of the 
position angles (13°) is significantly smaller than the half- 
opening angle of the X-ray caps, and the time sampling is 
currently too sparse to determine a precession period. Mul- 
tiple VLBI measurements of the position angle of the jets 
would be required to test the precession scenario. 



5 CONCLUSIONS 

We have resolved the milliarcsecond-scale jets in Circinus 
X-l for the first time. The jet morphology is symmetric, ex- 
tended on a scale of ~ 20 mas along a position angle 112° 
east of north, in good agreement with the position angle 
seen in ATCA observations of the arcsecond-scale jets. We 
do not resolve the jets along the minor axis, implying an 
opening angle of < 20°. The source brightness decayed over 
the course of the observing run, and there is evidence for 
outward motion of the components between the first and sec- 
ond halves of the observation, implying a separation speed 
of ~ 35 mas d -1 . This is greater than the proper motion es- 
timate derived assuming that the jets were launched at the 
onset of the radio flare at orbital phase zero, implying a pos- 
sible delay in the launching of the jets followi ng periastron. 
While the ultra-relativistic flow inferred by iFender et al.l 
|2004l ) cannot be definitively ruled out, the symmetry of 
the observed jets is difficult to reconcile with this scenario. 
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